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Abstract 



Using recently developed methods, we perform high-accuracy calculations 
of the susceptibility near (3 C for the D = 3 version of Dyson's hierarchical 
model. Using linear fits, we estimate the leading (7) and subleading (A) 
exponents. Independent estimates are obtained by calculating the first two 
eigenvalues of the linearized renormalization group transformation. We found 
7 = 1.29914073 ± 10~ 8 and, A = 0.4259469 ± 10~ 7 independently of the 
choice of local integration measure (Ising or Landau-Ginzburg). After a suit- 
able rescaling, the approximate fixed points for a large class of local measure 
coincide accurately with a fixed point constructed by Koch and Wittwer. 
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Scalar field theory has many important applications such as superconductivity, low en- 
ergy descriptions of quark-antiquark bound states and possibly the mechanism responsible 
for the mass generation of all the experimentally observed particles. However, there exists 
no approximate treatment of this theory which could pretend to compete in accuracy with 
quantum electrodynamics at low energy, where renormalized perturbation theory can be 
used to calculate the magnetic moment of the electron and the muon with more than ten 
significant digits. Given the inherent difficulties associated with the experimental probing 
of very short distances, accurate calculations compared with accurate "low-energy" experi- 
ments may become more of a standard procedure in the next century. 

In order to obtain a completely satisfactory treatment of scalar field theory in various 
dimensions, one needs an approximation scheme such that: a) the zeroth-order approxima- 
tion preserves the main qualitative features of the model, b) the zeroth-order approximation 
allows very accurate calculations, and c) the zeroth-order approximation can be improved 
systematically and in a way which preserves its computational advantages. We advocate 
here that hierarchical approximations such as the approximate recursion formula derived by 
K. Wilson |I[ or the related recursion formula which holds for Dyson's hierarchical model 
could be used as such a zeroth-order approximation. The fact that the approximate 
recursion formula satisfies the requirement a) is discussed at length in Ref. M. The fact 
that requirement b) is satisfied in the symmetric phase is explained in Refs. |3],|| and briefly 
reviewed below. The use of the hierarchical approximation solves some important problems 
encountered in practical calculations in lattice field theory. First, it possible to perform all 
the integrals appearing in the calculation of the zero-momentum Green's function in a much 
more efficient way than with the Monte-Carlo procedure. Second, the computing time scales 
only like the log of the number of sites, and one can eliminate finite-size effects completely. 
This can justify the effort of trying to solve part c) of the above program which is a very 
difficult problem. 

Wilson's approximate recursion formula is closely related to the recursion formula appear- 
ing in Dyson's hierarchical model ||. Both models have no wave function renormalization 
(77 = 0). It is possible to continuously interpolate between the two models and to show 
that during this process, the critical exponent associated with the susceptibility 7 varies || 
by less than 5 percent with respect to the nearest neighbor value. However, the numerical 
treatment of the two models is completely identical. In the following, we specialize the dis- 
cussion to the case of Dyson's model because this model has been studied [p|,|6|-|T0H in great 
detail in the past. We want to make clear that this choice is not essential and anything done 
below could have been done for Wilson's approximate formula. 

In a typical lattice field calculation, we pick some values for the bare parameters entering 
in an action and we calculate the renormalized quantities. In general, the physical masses 
cannot be too large when expressed in cut-off units. Ideally, we should be able to cover a 
broad range of situations going from effective theories with a low cut-off (e.g., (m p /m n ) ~ 6 
for an effective theory of pions) to a "fundamental" theory with a large cutoff and which 
requires some fine-tuning procedure. In the following calculations, the fixed bare parameters 
will appear in a local measure of the Landau- Ginzburg (LG) form: 

W ((f>) oc exp-^ m2<t,2+9 ^ ] . (1) 
The limit |l] of a large UV cut-off A, in units of the physical mass tur, can be reached by 
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tuning another parameter /3, which is the inverse temperature in Dyson's formulation of the 
model f2|. More explicitly, A/m.R is proportional to (/3 C — /3)~i . 

By fine-tuning f3, one can approach a fixed point of the renormalization group (RG) 
transformation and describe the RG flows using the linear approximation. As we will show 
later, for (3 close enough to (3 C (i.e. , for A large enough), one can approximate very well the 
magnetic susceptibility (zero-momentum two point function) with a linearized expression 
taking into account only the first irrelevant direction: 

x ^^ c -^)-\A Q + A l ^ c -^) , (2) 

If we were sure that there is only one non-trivial fixed point (universality) and if we we 
could calculate accurately the exponents, then the complicated procedure described above 
can be reduced to the determination of Aq and A\ in Eq.(|2]), a procedure that involves no 
fine-tuning. 

In this Letter, we provide empirical evidence suggesting that the RG transformation of 
Dyson's hierarchical model has only one non-trivial fixed point. We calculate the exponents 
7 and A with two independent methods (direct fit and linearization). The accuracy of 
our results is significantly better than the accuracy reached in the past All the 

approximate fixed points we have constructed below are very close (after rescalings explained 
below) to the fixed point calculated with an extraordinary accuracy by Koch and Wittwer 
||. Our work demonstrates the enormous calculational advantage of using the hierarchical 
approximation and addresses the question of understanding to what extent expansions about 
a known fixed point can be used as a substitute to the lengthy calculations in terms of bare 
parameters described above. 

For the sake of completeness, we briefly review the steps which lead to the basic expression 
of the RG transformation of Eq.(||). The block-spin transformation of the hierarchical model 
is an integral formula which transforms the local measure W(4>) according to the rule: 

W n+ M) <x e«<f" + '*' / W^V„(^±«) , (3) 

where c = 2 1- ^ in order to approximate D-dimensional nearest neighbor models. For more 
details, the reader may consult Refs. ||||||. In the following we only consider the case D = 3 
in the symmetric phase. We approach criticality for a fixed initial Wq by fine-tuning /3 as 
described in Ref. [|J. When a critical value [3 C is reached approximately, a (discrete) scale 
invariance is temporarily restored and it is convenient to reabsorb the scale factor (c/4) in 
<fr 2 . After this rescaling, we obtain the conventional RG transformation of the local measure. 
In Fourier form it reads: 

R n+ i(k) = C n+l eM-\^){Rn{^)Y ■ (4) 

We fix the normalization constant C n in such way that R n {0) = 1- We consider the finite 
dimensional approximations of degree l ma x' 

R n (k) = 1 + a n ^k 2 + a n , 2 /c 4 + + a nUax k 2lm ^ . (5) 

The coefficients a n j are proportional ||,[I0]|| to the expectation value of the sum of all the 
fields (after n iterations, there are 2 n of them) denoted M n . In particular, the finite volume 
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susceptibility Xn(P), denned as < (M n ) 2 > n /2 n is simply — 2a ri: i(^) n . When (3 < (3 C) Xn 
reaches a finite limit x when n goes to infinity. The recursion formula for a n m is purely 
algebraic: 



The initial condition for the Ising measure is Rq = cos(k). For the Landau- Ginsburg mea- 
sure, the coefficients in the k— expansion need to be evaluated numerically. 

In a recent article M ,we have shown that the errors on x due to finite volume and finite 
truncations fell exponentially fast with, respectively, the number of iterations used and the 
dimension of the truncated space (l ma x)- It is possible to make calculations where these 
errors play no practical role. The main limitation of the method comes from the round- 
off errors which are amplified when many iterations are spent near the fixed point. If the 
arithmetic operations are performed with a precision 5, then Q 

& 5 (7) 

l X l Pc-P' K> 

We now proceed to determine the values of the four parameters appearing in Eq. (|2|) from 
direct calculations of x a t various temperatures. The calculations which follow have been 
performed for two particular choices of Wo, one corresponding to the Ising limit (Wo(0) = 
5(<f) 2 — 1)) and the other to the choice m? — 1, p — 2, and g = 0.1 in Eq.(|l|). Unless 
specified differently, the calculations are performed using double-precision. In the following, 
we use the notation x for the quantity —logio(f3 c — j3). If we display log(x) versus x, the 
deviations from the linear behavior are not visible to the eye and need to be studied and 
understood "locally" in (3. In order to get a rough understanding of the corrections, we have 
divided the computer data in 14 bins of 100 points. The first bin contains data for values 
of x = 1.00, 1.01, . . . 1.99 and so on. In each bin (indexed i), we make a linear fit of log 10 (x) 
versus x. In the z-th bin, we call the slope 7W, and (cr^) 2 denotes the sum of the squares 
of the difference between the data and the linear fit divided by the number of points in a 
bin (100) minus 2. The values of a® are displayed in Fig. |1|. This graph can be interpreted 
easily. There are two known sources of deviations from the linear behavior: the subleading 
corrections to the scaling laws (which decrease when (3 gets close to (3 C ) and the round-off 
errors (which increase when (3 gets close to (3 C according to Eq.([?D). The approximate slopes 
in Fig. [I] confirm this interpretation. In bin 9, we minimize the combined deviations from 
linearity and we can consider 7^ as a first estimate of 7. Its numerical value is 1.29917 
in the Ising case and 1.29914 in the LG case. With this simple-minded procedure, we have 
already gained almost two significant digits compared to the existing estimates [0,0,0 where 
the answer 7 = 1.300 was consistently obtained with errors of order 1 in the last digit. 

We can improve this result by estimating the subleading corrections. For this purpose, 
we have used the bins 6 and 7 where the next subleading corrections are small (see discussion 
later) and where the numerical errors are still not too large. We have divided these two bins 
into 10 sub-bins of 100 points each. We use two digit indices for these sub-bins. For instance 
sub-bin 6.3 is the third sub-bin of bin 6 and contains the values of x: 6.3, 6.301, . . . ,6.399. 
Using Eq.(0), the same kind of notations as above for 7 and noting that j + 0.0495 is the 
middle of the sub-bin indexed by j, we obtain the approximate decay law: 
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7 a)^ 7 _ A( ^i)i ^a + o.o495) (8) 

The unknown coefficients can be extracted from linear fits of Zo(7io(7^ +0 — 7^'). We 
obtained 4 1 = -0.57 and A = 0.428 for the Ising model and 4 1 = 0.14 and A = 0.427 
for the LG model specified above. Repeating the first step (a linear fit in bin 9) but with 
X divided by (1 + j^((3 c — f3) A ), we obtain 7=1.299141 with an agreement up to the sixth 
decimal place between the two models considered above. 

Eq. (0) is an unavoidable limitation if we use double precision arithmetic. However, using 
Mathematica with a suitably set precision, l max = 42 for the Ising model and l max = 50 
for the LG model (see ref. Q for the determination of these quantities), we were able to 
calculate x i n Drns 10,11 and 12 with 11 significant digits. In the following, we call this 
data the "high-precision data" . Since this procedure is relatively lengthy, we have used only 
ten points per bins. We also determined (3 C with 24 significant digits so that in bin 12, the 
subtracted quantity j3 c — (3 is also known with at least 11 significant digits. In the Ising case, 
the result is easily reproducible and reads (3 C = 1.17903017044626973251189. We have then 
used bin 12 (where the sub leading corrections are very small and our errors on them are less 
important) with x divided by the subleading correction as described above, to estimate 7. 
We then used this better value of 7 to obtain the subleading corrections in bin 7 (where they 
are more sizable). This procedure can be iterated. This "bootstrap" of linear fits converges 
rapidly. We reach a nine significant digit agreement between the high-precision data and the 
fit obtained with the above procedure. The small discrepancies can be analyzed in terms of 
first order errors made in the estimate of the four parameters. This linear analysis provides 
small corrections (< 4 x 10~ 9 ) to 7 and more sizable corrections (< 4 x 10 -4 ) to A. The size 
of these corrections provide an order magnitude estimate for the errors. After these small 
corrections are taken into account, we obtain an agreement between the exponents of the 
two models for the following digits: 7 = 1.299140730 and A = 0.4260. We conclude that 
7 = 1.29914073 with an estimated error of less than 10~ 8 . 

We would like to comment about the corrections to Eq. (|2[) and how they could affect our 
estimates. First, since the third eigenvalue of the linearized RG transformation A3 ~ 0.48, 
the next subleading exponent is approximately 2. For x > 10, these effects are negligible. 
Second, a general argument 0, suggests that we should replace the constant A and A\ in 
Eq.(g) by log-periodic function which can be expressed as linear superposition of Fourier 

modes of the form (f3 c — P) lnl - X ^, with I an integer. Evidence for non-zero Fourier modes 
were found in Ref. M by using an estimator of 7 — 1 called the extrapolated slope and 
denoted S m . In this estimator, oscillating and constant contributions have roughly the 
same amplitude. However, using Eqs.(3.7) to (3.10) of Ref. [f|, one realizes that in S m , the 
oscillating amplitude is dramatically amplified by a factor of the order |a; 3 /r(7-|-2c<j)| where 
uj = 27r/Zn(Ai) ~ 18. This implies that the Fourier coefficients of the non-zero modes are 
suppressed by at least 14 orders of magnitude. A direct search for these oscillations confirms 
this upper bound. Third, Eq. (|2]) is obtained from a linearization. Higher order corrections 
give contributions proportional to (/3 C — /3) 2A . An analysis of the difference between fit and 
data in low bins indicates that these corrections are the main source of errors in our analysis. 

An alternative calculation of the exponents consists in linearizing the RG transformation 
near a fixed point. An approximate fixed point can be found by approaching f} c from below 
with our best resolution and iterating until a n +i,i/ a n,i takes a value which is as close as 
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possible to 1. In the present formulation, the linearized RG transformation is given by the 

Imax X l-max matrix 

M hm = (9) 

evaluated at the (approximated) fixed point. Using the high-precision Mathematica-b&sed 
method described above, we obtained this approximate fixed point for n = 101 for the Ising 
model and for n = 97 for the LG model. Calculating the eigenvalues of Eq.(||) for the 
two models used for the first estimates, we obtain discrepancies of 2 x 10~ 8 for Ai and of 
4 x 1(T 8 for A 2 . The average values are Ai = 1.42717246 and A 2 = 0.85941163. Changing 
n by one or improving the fixed point using Newton's method produce variations in these 
eigenvalues which are smaller than 3 x 10~ 8 . Using the relations 7 = 2ln(2)/3ln(Xi) and 
A = -ln(X 2 )/ln(X 1 ), we obtain 7 = 1.29914078 and A = 0.4259469 both with estimated 
errors of order 10~ 7 . The new estimate of 7 is compatible with the previous one but is less 
accurate. On the other hand, the new estimate of A is more accurate. The discrepancy with 
the previous estimate is less than 10~ 4 which is compatible with our previous error estimate. 

The two approximated fixed points obtained in the above calculation depend on (3 C . We 
denote them R*(k, (3 C ). However, it is possible to obtain what will turn out to be a universal 
function U (k) by absorbing (3 into k. More explicitly, we found that 

u(k) =iF(v&,&) (10) 

is in very good approximation independent of the model considered. This function is related 
to a fixed point f(s 2 ) constructed in Ref. || by the relation 

U(k) oc /((^V) • (11) 



The Taylor coefficients of / can be found in the file approx.t in ||. Normalizing Eq.flTTD 
with U(0) = 1, we obtain 

U(k) = 1. - 0.358711349882A; 2 + 0.053537288227& 4 - . . . (12) 

It is not known if there is only one non-trivial fixed point for Dyson's model. Both the two 
approximate fixed points we have constructed above give a function U (k) very close to Eq. 



(|12|). The closeness can be characterized by the p- norms introduced in ||. For p = 2 and 
Z < 42 we found that the error 5ui on the Z-th coefficients of the approximate U(k) with 
respect to the accurate expression obtained from Ref. || were bounded by \8ui\ < — . 

In order to further explore the possibility of having different fixed points, we have con- 
sidered more LG models. Using the parametrization of Eq. ([!]), we have considered the 12 
cases obtained by choosing among the following possibilities: m 2 = ±1 (single or double-well 
potentials), p = 2, 3 or 4 (coupling constants of positive, zero and negative dimensions when 
the cut-off is restored) and g = 10 or 0.1 (moderately large and small couplings). These 
searchs have been performed using regular double-precision calculations. We have not aimed 
at great accuracy. For all these twelve models, we found that using the same notations and 
conventions as a above \8ui\ < — . In other words, the function U{k) seems to be inde- 
pendent of the general shape of the potential, the strength of the interactions and whether 
or not the model is perturbatively renormalizable. 
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In conclusion, our best estimates of the critical exponents 7 = 1.29914073 ± 10~ 8 and, 
A = 0.4259469 ± 10~ 7 have an accuracy significantly better than existing estimates [|],[7|,[|. 
Our results demonstrate the power of the methods developed in Ref. M. They provide an 
incentive to develop more efficient perturbative calculations of the critical exponents and 
to attack the problem of the improvement of the hierarchical approximation. We found no 
indications for the existence of a non-trivial fixed point different from the one obtainable 
from Ref. M. Near criticality, or in field theoretical language for a large UV cut-off, the 
parametrization of Eq.(Q) fits the data very well. The quantities A and Ai depend on the 
bare parameters in a complicated way. However, the fact that we can use confidently the 
universal features suggests that it is possible to shortcut the use of bare parameters and 
consider directly Aq, A\ as an input. More generally, we are in position to check if the 
following conjecture is true: an expansion about the non-trivial fixed point can be used as a 
substitute for the calculations in terms of bare parameters. If true, this would mean that the 
result of Ref. || effectively "solves" the hierarchical model even far away from criticality. 

This research was supported in part by the Department of Energy under Contract No. 
FG02-91ER40664. J.J. Godina is supported by a fellowship from CONACYT. Y.M. thanks 
R Wittwer and the CERN lattice group for useful conversations. 
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FIG. 1. The deviations from the linear fits defined in the text as functions of the bins, for 
the Ising model (circles) and the Landau-Ginzburg model (stars). 
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